On the exact calculation of travelling wave solutions to nonlinear evolution equations 
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We analyse and compare three methods for the exact calculation of travelling wave solutions to 
nonlinear partial differential equations. We simplify the so called (G'/G)-expansion method and 
apply two of those methods to simple physical problems. 
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. There has recently been great interest in the application of the (G"/G)-expansion method to obtaining travelling 

. wave solutions of some nonlinear partial differential equations The method was earlier applied to a variety of 

' ^ '■ such problems^, There are many mathematical recipes for that purpose[ll, two of them are the exp-function 

JpH' method [5i and the tanh-function method [6]. All those prescriptions look quite similar and have also been applied to 

' similar toy problems, ft is curious that there has not been much interest in comparing them, 

g ' In this paper we analyze those recipes. In SeclUwe outline the (G"/G)-expansion method and in Sec. IIIII propose 

I— 1^ a simplified version of it. In Sec. IIVI we compare the three prescriptions just mentioned and in Sec. |V]we apply them 
to two simple examples. Finally, in Sec. IVII we draw conclusions from the results of the preceding sections. 

II. THE (G7G)-EXPANSION METHOD 
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\ This method was proposed to obtain analytical travelling-wave solutions to nonlinear differential equations of the 

0\ ; formQ 
O 



P{u,Ut,Utt,Uj;,Ua:x,Utx, ■ ■ ■) ^ (1) 

where the subscripts indicate differentiation of u{x,t) with respect to its arguments, and P is a polynomial function. 
The (G"/G')-expansion method has been applied to several examples of travelling waves that are particular solutions 
of the form u{x,t) — u(^), where ^ = a; — ct so that ut — ~cu' and = w' [l[ 3, [^j | . 

One assumes that the travelling-wave solution to the differential equation ([T]) is of the form 



40 = 



G(0 



(2) 



where the function G{^) is a solution to 

G"(O + AG'(O + AiG(O=0 (3) 



It follows from 
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2 



that d^u/d£^^ is a polynomial function of G jG of degree m + k. For this reason, substitution of Eq. ([2|) into Eq. ([T]) 
yields a polynomial function of G' /G and if we set its coefficients equal to zero we may obtain the unknown coefficients 
aj. 

III. THE (F7F)-EXPANSION METHOD 

The function F{(,) = e^^/'^G{0 is a solution of the differential equation F" + (/i — A^/4)i^ = and satisfies 

----- (5) 
G^ F 2 



Therefore the series ([2]) becomes 



F{0 



(6) 



where the coefficients bj absorb the parameter A that therefore dissapears from the equations. We may call 7 = /i— 
and simply require that the function F{^) is a solution to 

F"{O+^m = (7) 
For completeness and comparison we show the solutions to this trivial equation: 

{ci cos (/7^) + C2 sin (^7^) , 7 > 
ciCOsh(V^0 +C2sinh(V^C), 7 < (8) 
ci + 7 = 

Notice that this expression leads to all the cases explicitly shown by Zayed and GepreelQ. Therefore, this method is 
entirely equivalent to the (G"/G)-expansion method but the equation that defines the main function is simpler. In 
Sec.|V]we show an application of the (i^'/F)-expansion method. 

It is worth noticing that the simplification proposed here also applies to the generalized (G"/G)-expansion methodQ 
because it is based on the same equation for G(^). 

IV. COMPARISON OF DIFFERENT METHODS 

The function w{^) — ^F' {0/ F{0 is a solution to the Riccati equation 

w' = 7 + (9) 
proposed by FanQ to expand the traveling- wave solutions in the form 

m 

3=0 

We appreciate that the extended tanh-function method is equivalent in principle to the (i^'/_F')-expansion method 
and thereby to the (G'/G)-expansion method. There is a difference, however, in that the latter two methods appear 
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to be more flexible because they have the two independent solutions of the generating equation built in the expansion 
function through the constants ci and C2. 

We can rewrite the solution to Eq. ([7]) in a different way 



^^(0 = cie"« + C2e-"S a - (11) 

so that the expansion ([6]), and consequently also ([2]), reduces to a ratio of two polynomial functions of e"^ of degree 
2m. We thus conclude that the (G"/G')-expansion method is a particular case of the Exp-function methodlsl that is 
based on a solution of the form 

. . EU-c °J exp(j>?) 

"('y) = y4 Z -—-,ri^kx + ujt (12) 

Lj=_p 0] exp(j77) 

Notice that = if a = and c = —uj/k. It is clear that we can rewrite the rational approximation (|12p in such a 
way that it only shows positive powers of e'' in the numerator and denominator. 

In the following section we apply the Exp-function method to a problem that apparently cannot be treated by 
means of the (G'/G) -expansion method. 



V. APPLICATIONS 

It is not the purpose of this paper to abound with arbitrary tailor-made toy model equations for the application 
of the methods just described. However, in what follows we consider two well-known exactly solvable problems that 
will enable us to compare the methods outlined above. 

Our firs example is Fisher's equation 

du d'^u _ , 

that was originally derived for the simulation of propagation of a gene in a population and also arises in heat and 
mass transfer, biology, and ecology. This equation is so popular that appears in most on-line encyclopaedias of 
mathematics For our present purposes it suffices to mention the application of the Exp-function method^. 
If we choose £, = x — ct (we may also choose £^ — x + ct) we obtain 

u" + cu' + u{l - u) (14) 

The leading terms of u" and and are of dgree m + 2 and 2rn, respectively, and therefore they do not cancel unless 
m = 2. If we substitute Eq. © into Eq. ^ we easily obtain &2 = 6, fei = -6c/5, bo = [25(87 + 1) - c^J/SO, 
7 — — c^/100. The constant term gives us the only values of c for which there are exact solutions: c = ±5/V6- For 
c = 5/V6 we obtain 

uiO = 72 (15) 

1 + Ge«/^ 



By means of the Exp-function method Zhou[9| identified four cases that produced four solutions named Uj{x,t), 
j — 1, 2, 3, 4. However, they are not essentially different as one can show that they are closely related by the symmetry 
of the problem: ui{—x,t) = U2{x,t), U3{—x,t) = U4{x,t), and ui{x,~t) becomes U3{x,t) after substitution of 4/6o 
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for bo in the latter. On the other hand, Eq. (|15p is identical to U2{x,t) if we simply substitute 2/bo for C. The other 
solutions follow from the indicated symmetry or by choosing the other root c — — 5/V6- 
Our second example is the Bratu-Gelfand equation 

u"ix) + Ae"(^) = 0, u'{0) = = (16) 



that appears in simple models for the stationary theory of the thermal explosion [10[. First, we have to convert this 
strongly nonlinear equation into a polynomial one, which we do by means of the transformation u{x) = —nliiv{x) 
that leads to Au^~" + n (y'^ — vv"^ — 0. In order to have the simplest equation we choose n = 2, so that 

2 iy'^ - vv") + A = 0, v\Q) = 0, v{l) = 1 (17) 

If we substitute 

i'(a;) =a_ie-"^ + ao + aie"^ (18) 

which is a particular case of Eq. , into Eq. p7|) we easily obtain 

ao = 0, a_i = ^ „ (19) 

The boundary condtions w'(0) = and v[l) = 1 give us two additional equations 

A 

aia - = 

8aia 



+ ^ - 1 (20) 



respectively, from which we obtain 



(e2" + 1)^ 



«i - (21) 

The first equation shows a well-known bifurcation diagram from which one obtains the critical condition of ignition 
ckc = 1.19967864 that is a root of dX/da = o[l^. In other words, explosion takes place in a self-ignition process for a 
plane-parallel vessel when A = A^ = A(ac) = 0.8784576797Q. 

Neither the (G'/G')-expansion methodP, 3, 3, nor the tanh-function method^] provide a function of the form 
that is already given by the more flexible Exp-function method^]. 



VI. CONCLUSIONS 



There seems to be an ever increasing number of methods for the solution of more or less trivial nonlinear problems; 
many of them are cited in the papers by Zayed and Gepreel[l| and Ganji and AbdoUahzadehQ. Our analysis of three 
such recipes was motivated by the interest of this journal in one of those methods. We have shown how to simplify 
the main equation of the (G"/G)-expansion method and have compared it with the tanh-method and Exp-function 
methods. Our simple and straightforward discussion shows that the first two ones are similar in essence and that the 
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third one is more general and flexible. The three of them lead to solutions in the form of quotients of polynomial 
functions of exponential ones. 

We have illustrated the simplification of the (G'/G')-expansion method, named (F'/F)-expansion method, by its 
application to the simple and widely known Fisher's equation^, Q]. Finally, we have verified the greater flexibility of 
the Exp-function method with the aid of the Bratu-Gelfand equation 10]. 

The papers of Zayed and Gepreelilj and Gan ji a nd Abdollahzadeh ^ do not exhibit the amazing errors of other 
articles that I discussed in the past [111 However, they are standard applications of a method 

developed earlier and widely applied in previous papers by several authors (see also the references cited therein). 
It is suprising that Journal of Mathematical Physics (supposedly devoted to original and more elaborate mathematical 
applications to physics) accepts such kind of contributions. As an experiment I submitted a comment to JMP (with 
the content shown above) and it was rejected on the grounds that it did not "contain sufhciently significant new 
results to warrant its publication in JMP" . Does it seem that the Elsevier policy of publishing poor papers for reasons 
other than purely scientific ones is beginning to contaminate also JMP? 
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